Texture-induced modulations of friction force: the fingerprint effect 
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Dry solid friction is often accompanied by force modulations originating from stick-slip instabili- 
ties. Here a distinct, quasi-static mechanism is evidenced leading to quasi-periodic force oscillations 
during sliding contact between an elastomer block, whose surface is patterned with parallel grooves, 
and finely abraded glass slides. The dominant oscillation frequency is set by the ratio between the 
sliding velocity and the period of the grooves. A mechanical model is proposed that provides a 
quantitative prediction for the amplitude of the force modulations as a function of the normal load, 
the period of the grooves and the roughness characteristics of the substrate. The model's main 
ingredient is the non-linearity of the friction law. Since such non-linearity is ubiquitous for soft 
solids, this "fingerprint effect" should be relevant to a large class of frictional configurations and 
might in particular have important consequences in human (or humanoid) active digital touch. 

PACS numbers: 46.55.-|-d, 68.35.Ct, 81.40.Pq 
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Surface texture engineering by nano- or micro- 
patterning techniques has proven to be an efficient tool to 
tune frictional or adhesiveproperties at solid/solid ll|-|3[ 
or solid/liquid interfaces [J|. Many of the existing de- 
signs arc directly inspired by natural surfaces. Super- 
hydrophobic surfaces for example, reproduce the micro- 
scopic pattern observed on lotus flower leaves [5|. Sur- 
faces with a shark skin-like texture provide significant 
drag reduction in turbulent flows [6[. Geckos feet fib- 
rillar structure provides them with remarkable adhesion 
capabilities, and has allowed the design of functional ad- 
hesives [7|. 

In comparison, the effect of micro-patterning in 
solid/solid friction has so far attracted much less atten- 
tion. In a recent paper however, we suggested that epi- 
dermal ridges - the regular pattern characteristics of hu- 
mans and primates glabrous skin surface - strongly mod- 
ify the skin internal stress elicited upon actively rubbing 
the fingertip over a finely abraded surface, thus enhanc- 
ing tactile sensitivity [8|. More precisely, such patterns 
were shown to produce a quasi-periodic modulation of 
the texture-induced subcutaneous stress field at a fre- 
quency set by the ratio between the scanning velocity 
and the pattern period. This finding led to the integra- 
tion of fin gerp rint-like structure in several robotic hand 
designs |9|, llO| . Interestingly, this effect not only shows 
up in local stress measurements but is also observed in 
the global friction force signal measured while scanning 
a human finger on a rough surface [llj . Since the friction 
force signal conveys a large part of the tactile informa- 
tion for texture discrimination, as recently evidenced in 
psychophysical assays |12| , it is crucial from both biologi- 
cal and robotic point of views to understand the physical 
origin of such texture-induced friction force modulations 
and to identify the parameters which control their spec- 
tral characteristics. 

To address this question, we investigate in this 
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FIG. 1. (color on line) Sketch of the experimental setup. 
Inset: image of a typical contact at Fn = 0.5 N with 
A = 125 ^m. The red circle defines the apparent contact 
perimeter. 



Letter the friction dynamics between a patterned elas- 
tomer block and a finely abraded glass surface in a 
sphere-on-plane geometry (Fig.[T]). The elastomer blocks 
(thickness «18 mm) were obtained by molding liquid 
PolyDiMethylSiloxane-crosslinker mixture (PDMS, Syl- 
gard 184, Dow Corning) in concave spherical lenses of 
radius of curvature R = 128.8 mm. Prior to molding, 
the lenses surface was patterned with a periodic square 
grating of depth 40 ^m and spatial period A = 125, 218 
or 760 /im (width A/2) using soft photolithography tech- 
niques [8|. The elastomer Young's modulus was mea- 
sured to be S = 2.2 ±0.1 MPa. The abraded planes con- 
sisted in microscope glass slides sandblasted or abraded 
using a Silicon Carbide powder-water mixture. The sur- 
face topography was measured with an optical profilome- 
ter (M3D, Fogale Nanotech) from which the root mean- 
squared height hrms and surface profile autocorrelation 
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FIG. 2. (color on line) (a)-(d) Blow up oiFsjFs fluctuations 
(resp. Fn /Fn in thin black lines) versus u/\ at Fjv=0.5 N 
for different A and both substrates in steady sliding, (e) Av- 
eraged power spectra of Fs fluctuations (A=125 /im, Fjv=0.5 
N, "rough-" substrate), both with v perpendicular (•) and 
parallel (o) to the ridges. 



spectrum C{q) were computed [13|. Two substrates were 
considered, referred to as "rough-" {hrms ~ 2.3 ^m) 
and "rough-f" (hrms ~ 5 ^im), respectively. For both, 
C{q) displays a self-afBne behavior as usually observed for 
real surfaces [1J| of the form C{q) = Cq (l + [q/qcY) 
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with (Co,gc,a) 
(1.3 10-20 ni4, 8 
rough+, respectively. Substrates were rubbed against the 
elastomer surface along the x direction at constant veloc- 
ity V (range 0.01 - 0.5 mm.s~^) with a motorized trans- 
lation stage (LTA-HL, Newport), under constant normal 
force Fm (range 0.02 - 2N). Fm and tangential force 
Fs were measured with a 1 mN accuracy at 1 kHz by 
monitoring the deflections of two orthogonal cantilevers 
(Fig. [1]) using capacitive position sensors (MCC-10 and 
MCC-20, Fogalc Nanotech). The apparent contact zone, 
imaged in a light-transmitted geometry (Fig. [T|) had a 
diameter ranging from 1 to 6 mm, much larger than the 
substrate roughness characteristic scale. 

Typical Fs signals in steady sliding are shown on 



Fig. [2^-d as a function of the scanned distance u = v.t 
with V = 0.1 mm.s~^ perpendicular to the ridges direc- 
tion and Fjv = 0.5 N. Minute but clearly measurable 
quasi-periodic oscillations are observed, with no equiva- 
lent in the normal force signals. Their amplitude weakly 
depends on A and increases with Fjy. They are signifi- 
cantly more pronounced with the rough+ than with the 
rough- substrate (Fig. ^p and [2Jl) . The effect entirely 
vanishes when the ridges direction is aligned with the di- 
rection of motion. In the range of velocities explored, 
these modulations are independent of v. All experiments 
presented further are thus done aX v ~ 0.1 mm.s^^ and 
u is taken as the time- varying variable. 

The power spectrum 5'[5i^g] (g) of 5Fs {u) = Fs {u) — Fs 
(the bar stands for time averaging) was computed by av- 
eraging over a 20 mm travel distance. All spectra exhibit 
well defined peaks at q\ ~ 27r/A and at corresponding 
harmonics (Fig. [2j;). S'l^i^gj a,t q = qx increases with _F)v 
as AF'^ (Fig. El where values of A and v are collected 
in Table 1. As mentioned above, A depends weakly on 
A and is significantly higher for the rough-|- than for the 
rough- substrate. The exponent v is close to 1 for all 
experimental configurations. 

In parallel to the force measurements, rapid imaging 
of the contact zone was performed using a fast camera 
(Fastcam APX-RS, Photron, Japan) operating at 60 
Hz. At all loads, contact occurs only at the ridges 
summits, thus yielding a large contrast between the top 
and bottom of the pattern and allowing for a precise 
tracking of the ridges edges with « 10 /im accuracy. No 
stick-slip motion of the ridges was observed. Thus, the 
measured force oscillations cannot be accounted for by 
periodic stick-slip events whose frequency is set by the 
period of the pattern as reported in [2| . This is further 
supported by the observation that the fluctuations are 
quasi-static (u independent). 



In the rest of this Letter, we aim at understanding the 
physical origin of the observed friction force modulation 
and the dependence of its amplitude with the normal 
load, pattern period and substrate roughness character- 
istics. We flrst consider a perfectly smooth elastomer 
surface (no patterning) rubbed against an abraded glass 
slide. We denote p{x,y,u) and T{x,y,u) the normal and 
tangential stress field at a given position u of the sub- 
strate, such that Fn{u) = J p(x,y,u) dxdy and Fs{u) = 
J t(x, y, u) dx dy. Under a local Amontons-Coulomb fric- 
tion law assumption i.e. t{x, y, u) = ii,op{x, y, u) with ^o 
a uniform (material dependent) friction coefficient, the 
friction force reads Fs{u) = [J-oFn and is thus expected to 
be time-independent. We therefore hypothesize a weakly 
non-linear relationship between p{x,y,u) and T{x,y,u) 
which is equivalent to postulating a pressure dependence 
of the local friction coefficient, i.e. ^ = r/p = ^i{p). The 
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FIG. 3. (color on line) (a)-(d) SisPslil^) versus Fn- Dashed 
lines are best experimental fits, solid lines are model predic- 
tions. Shaded areas bound the values of A by tuning the 
parameters of the model. 



friction force now reads 

Fs{u) = / fJ.{p)pix,y,u)dxdy 



(1) 



The spatial variations of p{x, y, u) can be decomposed 
into a time-averaged component p{x, y) set by the macro- 
scale plane-on-sphere contact geometry and a time- 
fluctuating component 5p{x, y, u) associated with the mi- 
croscopic roughness of the glass slide. For simplicity, we 
further assume that the elastomer is in intimate contact 
with the glass surface. This assumption is expected to 
fail at the micro-asperities scale but should be valid at 
intermediate scale such as the inter-ridge distance A. In 
this limit, the texture-induced pressure modulations are 
set by the topography of the glass substrate and are in- 
dependent of the local mean pressure imposed by the 
Hertzian geometry p{x, y) . The pressure modulation field 
Sp(x,y,u) is thus a sole function of the position of the 
substrate with respect to the contact zone (u — x,y), so 
that the pressure field reads 

p{x, y, u) = K(u) {p{x, y) + 5p{u - x, y)) (2) 

The normalization factor K{u) ensures that Fj^ remains 
constant during sliding. This condition imposes 

K{u) = f 1 + — / 5p{u - x, y) dx dy 



1 - — / Sp{u - X, y) dx dy 



(3) 



TABLE I. : Experimental and model parameters. 
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limit, the integral of the roughness-induced pressure field 
becomes vanishingly small with respect to the confining 
force Fn- Rewriting Eq. [1] one obtains 

5Fs{u) = Fs{u) - Fs= I {fi{p) - {i^L))Sp{u - X, y) dxdy 

where (/i) is defined as the ratio Fs/Fn- 

The friction force fiuctuations thus appear as the con- 
volution product of a function characterizing the fric- 
tion coefficient spatial heterogeneities and the texture- 
induced pressure modulations field. The presence of reg- 
ular ridges at the surface of the elastomer imposes that 
the stress between ridges vanishes. This can be accounted 
for, in first approximation, by introducing a A-periodic 
Heavysidc function H{x) = 0{sin{2iTx / X)) under the in- 
tegral, which directly results in a spectral selection of the 
associated spatial mode. 

In Fourier space [l3| , the fiuctuating component of the 
friction force thus reads 



S^SFs]{qu) = (27r)4 / dq, \F{H{x){ii{x,y) - (m)))|'S| 



The latter expansion is valid when the contact diameter 
is much larger than the texture scale. Indeed, in this 
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Wave vector q^ (resp. g.„) refers to the parallel (resp. 
perpendicular) direction to the sliding direction. 

Under the assumption of intimate contact, Sr^pi can 
be simply expressed as a function of the roughness spec- 
trum C{q) [ll as S[sp\{q) = {E* /2f q^C{q), where q is 
the norm of {qu,qv) and E* the reduced Young's mod- 
ulus. The friction coefficient spatial field, on the other 
hand, is estimated from global force measurements. For 
all experimental configurations, a power-law dependence 
is observed between both force components Fs = B.F^ 
with 7 = 0.87 ± 0.04 (Fig. |4]). We postulate a power-law 
relationship between the local shear stress and the local 
pressure, T{x,y) = /3p^{x,y) [l5|. Assuming that p{x,y) 
follows a Hertz profile, the friction force Fs can be de- 
rived by analytically integrating the latter equation over 
the contact area. This procedure yields approximated 
values for /3 and m as function of B and 7. Figure |4] shows 
a comparison between the Fs versus F^ obtained exper- 
imentally and through integration of t{x, y) = (3p"^(x, y) 
over the contact area. 

Both predicted and measured amplitudes of S[sp^-^{q\) 
with Fn are presented on Fig. [3l for all experimental 
configurations, using the set of parameters (/3, m, Co, qc, 
a) obtained as described earlier. The predicted ampli- 
tude of Fs modulations is in reasonable agreement with 
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FIG. 4. (color on line) Time- averaged friction force Fs versus 
Fn- Different symbols correspond to different configurations. 
Inset: Average dynamic friction coefficient Fs/Fn versus Fn- 
On both graphs, lines are power law fits with 7 — 0.87 ±0.04. 
The dashed line has a slope of 1. 



the data, without any adjustable parameter. The model 
captures in particular the analytic increase of S[sFs]{<1\) 
with Ff^, with an exponent v comparable to its experi- 
mental counterpart (see Table 1 and Fig. [3]). By tuning 
systematically the input of the model (to, A), i^ is found 
to increase linearly with m and to decrease linearly with 
A in the investigated range (not shown). The large un- 
certainty of ly (~ 0.1) in the model is thus mostly due to 
the experimental inaccuracy in evaluating the parameter 
m. The predicted amplitude A is trivially proportional 
to Co, slightly decreases with A and falls rapidly to zero 
when 771 approaches unity, i.e. when Amontons-Coulomb 
friction's law is verified locally. 

This study shows that any non-linearity in the fric- 
tion law leads to the development of texture-induced 
friction force fluctuations. Since purely linear friction 
laws are scarcely observed, the present mechanism should 
be relevant to most practical situations. However, such 
fluctuations are expected to be in practice hardly de- 
tectable when the fixed substrate (here the PDMS block) 
is smooth and the contact zone diameter is much larger 
than the typical roughness scale. In contrast, the pres- 
ence of a regular pattern at the surface of the block op- 
erates a spectral filtering of the texture induced pressure 
modulation which induces a mode selection of the force 
fluctuations at a particular frequency. The amplitude 
of these force fluctuations can then become an order of 
magnitude larger than in the smooth case as shown by 
the present study. This mechanism was here examined 
for the simplest possible pattern, i.e. regular parallel 
stripes, but richer force signal spectra should be simi- 
larly obtained with more complex patterns. 

In essence, the macroscopic friction force signal car- 
ries information about the spectral content of the micro- 
scopic substrate topography at frequencies defined by the 
elastomer micro-pattern. A practical consequence is that 
Fs fluctuations amplitudes could be used to discriminate 



surfaces having small differences in their roughness at 
scales much smaller than the contact extension. Hence, 
under the same conditions, force fluctuations are roughly 
10 times larger for the substrate having a hrms ~ 5^m 
than for the one with hrms ~ 2.3/xm. This mechanism 
could also facilitate the detection of motion. The transi- 
tion between static and sliding contact should manifest as 
the sudden appearance of large single mode oscillations 
of the friction force. One may finally foresee the pos- 
sibility to evaluate the scanning velocity, based on the 
sole friction force signal, by continuously extracting the 
peak frequency of these oscillations. We believe that this 
"spcctrotribometry" approach may be relevant to human 
digital touch and prehension, owing to the presence of 
fingerprints, and that it could be easily implemented in 
tactile robotic sensing devices. 

The authors deeply thank F. Petrellis for illuminat- 
ing mathematical surgery, F. Zalamea for his valuable 
help, A. Chateauminois and C. Fretigny for fruitful dis- 
cussions, and acknowledge financial support from ANR- 
DYNALO NT09-499845. 



[1] 
[2] 

[s; 
[4; 
[5; 
[6; 

[7 

[s; 
[9; 

[10 

[11 

[12 

[13 

[14 
[15 



|alexis.prevost@lps.ens.fr| 

R. J. Charles and A. J. Crosby, J. Appl. Phys. 106, 

064913 (2009). 

R. Bennewitz, J. David, C.-F. de Lannoy, B. Drevniok, 

P. Hubbard-Davis, T. Miura, and O. Trichtchenko, J. 

Phys.: Conders. Matter 20, 015004 (2008). 

F. Wu-Bavouzet, J. Cayer-Barrioz, A. Le Bot, 

F. Brochard-Wyart, and A. Buguin, Phys. Rev. E 82, 
031806 (2010). ' 

P. Joseph, C. Cottin-Bizonne, J. M. Benoit, C. Ybert, 

C. Journet, P. Tabeling, and L. Bocquet, Phys. Rev. 
Lett. 97, 156104 (2006). 

B. Bhushan, Y. C. Jung, and K. Koch, Phys. Trans. 

Royal Society A-Math. Phys. and Eng. Sciences 367, 

1631 (2009). 

B. Dean and B. Bhushan, Phil. Trans. R. Soc. A 368, 

4775 (2010). 

L. F. Boesel, C. Greiner, E. Artz, and A. del Campo, 

Adv. Mat. 22, 21252137 (2010). 

J. Scheibert, S. Leurent, A. Prevost, and G. Debregeas, 

Science 94, 262 (2009). 

Y. Zhang, J. Micromech. Eng. 20, 085012 (2010). 

D. D. Damian, H. Martinez, K. Dermitzakis, 
A. Hernandez-Arieta, and R. Pfeifer, lEEE/RSJ 
2010 Int. Conf. Int. Robots & Syst. (IEEE, 2010) pp. 
904-909. 

A. Prevost, J. Scheibert, and G. Debregeas, Commu- 
nicative and Integrative Biology 2, 422 (2009). 

M. Wiertlewski, J. Lozada, and V. Hayward, IEEE 
Trans. Robotics 27, 461 (2011). 

B. N. J. Persson, J. Phys.: Conders. Matter 20, 312001 
(2008). 

B. N. J. Persson, J. Ghem. Phys 115, 3840 (2001). 

G. H. Scholz, The Mechanics of Earthquakes and Faulting 
(Cambridge University Press, 2010). 



